Evidence for the charmless annihilation decay mode 5° — t- tt^ 
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We search for annihilation decay modes of neutral b mesons into pairs of charmless charged hadrons 
with the upgraded Collider Detector at the Fermilab Tevatron. Using a data sample corresponding 
to 6 fb~^ of integrated luminosity, we obtain the first evidence for the — >■ ■k'^-k~ decay, with 
a significance of 3.7(t, and a measured branching ratio B{B'l — >■ -k^-k') = (0.57 ± 0.15 (stat) ± 
0.10 (syst)) X 10"®. A search for the B° -5> K+ K' mode in the same sample yields a significance 
of 2.0(7, and a central value estimate B{B° K+K') = (0.23 ± 0.10 (stat) ± 0.10 (syst)) x 10"®. 



PACS numbers: 13.25.Hw 14.40. Nd 

Our understanding of the dynamics of hadrons con- 
taining heavy quarks has made great progress in recent 
years. The development of effective theories has allowed 
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increasingly accurate predictions for the partial decay 
widths of such hadrons. An ability to make accurate 
predictions for these processes is not only important in 
itself, but is a tool to uncover possible additional contri- 
butions due to interactions beyond the standard model. 
In spite of the general progress of the field, a specific 
class of decay amplitudes (annihilation topologies) has 
resisted attempts at quantitative prediction up to the 
present, and is often simply neglected in calculations. 
Predictions for these amplitudes vary greatly between 
approaches, and even within the same approach. Es- 
timates based on the QCD factorization (QCDF) ap- 
proach are affected by significant uncertainties, due to 
end-point singularities More recent perturbative 

QCD calculations (pQCD) provide more precise predic- 
tions, but they tend to be significantly larger than the 
predictions coming from QCDF Q [if. No calculations 
are yet available within the soft collinear effective the- 
ory (SCET) [5|]. The lack of knowledge of the size of 
annihilation-type amplitudes introduces irreducible un- 
certainties in the predictions for several decays of great 
interest in the search for new physics effects, such as 
B° TT+TT- and B° K+K' (B-Q. Experimental 
investigation of the issue is therefore very desirable, and 
has the potential to enable a significant advancement of 
the field. The B° tt+tt" and B° K+R- decay 
modes are ideal for this investigation, because all quarks 
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FIG. 1: PA (left panel) and E (right panel) diagrams con- 
tributing to B° K*K~ and B° tt+tt" decays. 



in the final state are different from those in the initial 
state, so they can be mediated solely by amplitudes with 
penguin-annihilation {PA) and PF-exchange {E) topolo- 
gies (see Fig. [T]). However, they have not yet been ob- 
served, the best upper limits at 90% CL being respec- 
tively 1.2x 10"^ [lal and 0.41x10-^ A simultaneous 
measurement of branching fractions of both modes would 
be especially useful, as it would allow a better constraint 
on the strength of PA and E amplitudes W\ ■ 

In this Letter we report the results of a simultane- 
ous search for the two decays tt^tt" and 5° — > 
K^K~ using data corresponding to 6 fb^^ inte- 
grated luminosity of pp collisions at ^/s = 1.96 TeV, col- 
lected by the upgraded Collider Detector (CDF II) at the 
Fermilab Tevatron. 



The CDF II detector is described in detail in Ref. [13 1 
with the detector subsystems relevant for this analysis 
discussed in Ref. 14| . The data are collected by a three- 
level online event-selection system (trigger). At level 1, 
tracks are reconstructed in the transverse plane [l5| . 
Two opposite-charge particles are required, with recon- 
structed transverse momenta pti,PT2 > 2 GeV/c, the 
scalar sum pxi + Pt2 > 5.5 GeV/c, and an azimuthal 
opening angle < 135°. At level 2, tracks are com- 
bined with silicon-tracking-detector hits and their im- 
pact parameter d (transverse distance of closest approach 
to the beam line) is determined with 45 /im resolu- 
tion (including the beam spread) and required to be 
0.1 < d < 1.0 mm. A tighter opening-angle require- 
ment, 20° < A0 < 135°, is also applied. Each track pair 
is then used to form a B candidate, which is required 
to have an impact parameter ds < 140 /im and to have 
traveled a distance Lt > 200 /im in the transverse plane. 
At level 3, a cluster of computers confirms the selection 
with a full event reconstruction. 

The offline selection is based on a more accurate de- 
termination of the same quantities used in the trigger, 
with the addition of two further observables: the isola- 
tion (Ib) of the B candidate [31, and the quality of the 
three-dimensional fit (x^ with 1° degree of freedom) of 
the decay vertex of the B candidate. Requiring isolated 
candidates further reduces the background from light- 
quark jets, and a low reduces the background from 
decays of different long-lived particles within the event, 
owing to the good resolution of the silicon-tracking de- 



tector in the z direction. We use the same final selection 
originally devised for the — )■ K~Tr~^ search |10j|, whose 
simulation has proven to be nearly optimal also for detec- 
tion of — 7r+7r~. This includes the following criteria: 
Ib > 0.525, < 5, d > 120 fj.m, dB < 60 fim, and 
Lt > 350 /im. 

At most one B candidate per event is found after this 
selection, and a mass (m^+^-) is assigned to each, using 
a charged pion mass assignment for both decay prod- 
ucts. The resulting mass distribution is shown in Fig. [21 
and is dominated by the overlapping contributions of 
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17|, with backgrounds coming from misre- 



the B° - 
modes [li 

constructed multibody 6-hadron decays (physics back- 
ground) and random pairs of charged particles (combi- 
natorial background). A B'^ — > K'^K~ signal would 
appear in this distribution as an enhancement around 
5.18 GeV/c^, while a S° — > tt+tt^ signal is expected at 
the nominal mass of 5.3663 GeV/c^, where other more 
abundant modes also contribute [101] . 

We used an extended unbinned likelihood fit, incorpo- 
rating kinematic (kin) and particle-identification (PID) 
information, to determine the fraction of each individual 
mode in the sample. The likelihood is defined as 
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where N is the total number of observed candidates, ly is 
the estimator of N to be determined by the fit, and the 
likelihood for the ith event is 
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where the index j runs over all signal modes, and the 
index 'p' ('c') labels the physics (combinatorial) back- 
ground terms. The fj are the signal fractions to be deter- 
mined by the fit, together with the background fraction 
parameters b and /p. 

For each charged hadron pair, the kinematic informa- 
tion is summarized by three loosely correlated observ- 
ables: the squared mass m^+^_ ; the charged momentum 
asymmetry /3 = (p+ — p_ ) / (p+ + p_ ) , where p+ (p- ) is 
the momentum of the positive (negative) particle; and 
the scalar sum of particle momenta ptot — P+ + P~ ■ The 
above variables allow evaluation of the squared invari- 
ant mass m^+jj_. of a candidate for any mass assignment 
of the positive and negative decay products (ma+,mf,-), 
using the equation 



m 



where p^ 



Ptot^ , P~ 



1-/3 



m 



■ml. 12,) 



5 



> 

03 



0) 
CL 
0) 
Q) 
« 

■g 

« 
O 




I B^^ KV 
I B?^ TZ*n 



Multibody B decays 
Combinatorial bl<g 



5.6 5.8 



FIG. 2: Mass distribution of reconstructed candidates. The 
charged pion mass is assigned to both tracks. The sum of the 
fitted distributions and the individual components of signal 
and background are overlaid on the data distribution. 



olution tails are accounted for by a full simulation of the 
detector, while the effects of soft photon radiation in the 
final state are simulated by PHOTOS 2l| . This resolution 
model was accurately checked against the observed shape 
of the 3.2 X 10^ D° K-tt+ and 140 x 10^ tt+tt" 
signals in a sample of _D*+ — > D'^tt^ decays, collected 
with a similar trigger selection. As a result, the sys- 
tematic uncertainty related to the signal mass shapes is 
negligible with respect to other uncertainties. 

The D*'^ — )• D^TT^ sample was also used to calibrate 
the dE I dx response of the drift chamber to kaons and 
pions, using the charge of the D*+ pion to identify the 
I?" decay products. The dE / dx response of protons was 
determined from a sample of about 167 000 A — >■ p-K~ de- 
cays, where the kinematic properties and the momentum 
threshold of the trigger allow unambiguous identification 
of the decay products PID information is summa- 

rized by a single observable k, defined as: 



dE/dx - dE/dx{TT) 
dE/dx{K) - dE/dx{Ti) 



(4) 



TABLE I: Yields and significances of rare mode signals. The 
first quoted uncertainty is statistical; the second is systematic. 



Mode 




Significance 


B" K+K- 


120 ± 49 ± 42 


2.0(7 




94 ± 28 ± 11 


3.7cr 



The likelihood terms describe the kinematic dis- 
tributions of m^+^_ , P, and ptot variables for the physics 
signals and are obtained from Monte Carlo simulations. 
The same distributions for the combinatorial background 
are instead extracted from real data 18], and are inserted 
into the likelihood through the £j^'" term. In particu- 
lar, the squared-mass distribution of the combinatorial 
background is parametrized by an exponential function. 
The slope is fixed in the fit to the value extracted from 
an enriched sample of two generic random tracks, con- 
taining events passing all requirements of final selections 
except for vertex quality, replaced by an antiselection 
cut > 40, which strongly rejects track pairs originat- 



The likelihood term C^" 



ing from a common vertex, 
describes the kinematic distributions of the background 
from partially reconstructed decays of generic B hadrons. 
The m^j. _ distribution is, in this case, modeled by an 
ARGUS function [l9| convoluted with a Gaussian reso- 
lution, while /? and ptot distributions are obtained from 
Monte Carlo simulation. 

The fit has 28 free parameters. A detailed description 



of the fit and its parameters can be found in Ref. [IE 

To ensure the reliability of the search for small signals 
in the vicinity of larger peaks, the shapes of the mass 
distributions assigned to each signal have been modeled 
in detail. Momentum dependence and non-Gaussian res- 



where dE/dx{'K) and dE/dx{K) are the expected dE/dx 
depositions for those particle assignments. The average 
values of k expected for pions and kaons are by construc- 
tion and 1. Statistical separation between kaons and 
pions is about 1.4cr, while the ionization rates of protons 
and kaons are quite similar in the momentum range of 
interest. The PID likelihood term, which is similar for 
physics signals and backgrounds, depends only on k and 
on its expectation value (k) (given a mass hypothesis) 
of the decay products. In particular the physics signals 
model is described by the likelihood term >Cj^^, where 
the index j uniquely identifies the final state, while the 
background model is described by the two terms Cp^^ 
and C^^^, respectively for the physics and combinatorial 
background, that account for all possible pairs that can 
be formed combining only pions and kaons. In fact muons 
are indistinguishable from pions with the available dE/dx 
resolution, and are therefore included within the nominal 
pion component. For similar reasons, the small proton 
component in the background has been included within 
the nominal kaon component. Thus the physics back- 
ground model allows for independent, charge- averaged 
contributions of pions and kaons, whose fractions are de- 
termined by the fit; while the combinatorial background 
model, instead, allows for more contributions, since in- 
dependent fractions of positively and negatively charged 
pions and kaons are determined by the fit. 

The signal fractions returned by the fit are in agree- 
ment with those obtained in the previous iteration of 
this analysis [13]. The yields for the B'^ — > tt+tt^ and 
B'^ — > K^K~ modes, obtained from those fractions, are 
shown in Table [H The significance is evaluated as the 
ratio of the yield observed in data to its total uncer- 
tainty (statistical and systematic uncertainties added in 
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FIG. 3: Distribution of the relative signal likelihood, 
£-s/{Cs + Cothcv), in the region 5.25 < m^+^- < 5.50 GeV/c^ 
for B° 



-'Other ) 

K+K~. 



and 5.10 < m^+^- < 5.35 GeV/c^ for 
For each event, Cs is the likelihood for 



the B° TT+TT" (top panel) and B° K+K~ (bottom 
panel) signal hypotheses, and jCotha is the likelihood for ev- 
erything but the chosen signal, i.e., the weighted combination 
of all other components according to their measured fractions. 
Points with error bars show the distributions of data and his- 
tograms show the distributions predicted from the measured 
fractions. Zoom of the region of interest is shown in the inset. 



quadrature), where the statistical uncertainty is deter- 
mined from a simulation where the size of that signal is 
set to zero. This evaluation assumes a Gaussian distribu- 
tion of yield estimates, supported by the results obtained 
from repeated fits to simulated samples. This procedure 
yields a more accurate measure of significance than the 
purely statistical estimate obtained from y/ — 2Aln(£). 

We obtain a 3.7(7 significant signal for the B'^ — tt+tt" 
mode, and we observe an excess at the 2. Oct level for the 
B'^ — > K^K^ mode. As a check on the method. Fig. [3] 
shows relative likelihood distributions for these modes, 
which are in good agreement with our model. 

As a further check an alternate fit was performed, using 
kinematic information only. Removal of dE / dx informa- 
tion leads to results in agreement with the main fit, but 
with a loss in resolution of a factor 2 for i?g — > tt+tt^ 
and 3 for B^ — > K^K^, confirming the importance of 
this information. 

To avoid large uncertainties associated with produc- 
tion cross sections and absolute reconstruction efficiency, 
we measure all branching fractions relative to the i?° — >■ 
K+t:- mode. A frequentist limit [13] at the 90% C.L. 
is quoted for the i3° — >■ K~^K~ mode. The raw frac- 
tions returned by the fit are corrected for the differences 
in selection efficiencies among different modes, which do 
not exceed 10%. These corrections are determined from 
detailed detector simulation, with only two exceptions 



that are measured from data: the momentum-averaged 
relative isolation efficiency between 5° and 5°, and the 
difference in efficiency for triggering on kaons and pions 
due to the different specific ionization in the drift cham- 
ber. The former is determined as 1.00 ± 0.03 from fully- 
reconstructed samples of ^ and B*'— > J/ipK*^ 
decays jl^. The latter is determined from samples of 
mesons decaying into pairs of charged hadrons [isf . 
We measure the relative branching fractions S(D° — >■ 
Ti+TT-)/B{D^ ^ K-TT+) and ^(1?° ^ K+K-)/B{D° 
K~7:^). The numbers of events are extracted from the 
available samples of tagged tt+tt^, K^tt^ 
and — K^K~ decays, fitting the invariant D*tt mass 
spectrum [18| . while reconstruction efficiencies are de- 
termined from the same simulation used for the mea- 
surements described in this Letter. Comparison of these 
numbers with world measurement averages (23| allows us 
to extract the correction needed to compensate for the 
different efficiency of the tracking trigger for kaons and 
pions. The final corrections applied to our result do not 
exceed 5% and are independent of particle momentum. 

The dominant contribution to the systematic uncer- 
tainty on both branching fractions is due to the dE / dx 
model, which derives from the statistical uncertainty on 
the 48 parameters used for the analytical description of 
the correlated dE / dx response of the two decay prod- 
ucts [2^. This uncertainty is evaluated by repeating the 
likelihood fit 200 times with different sets of those pa- 
rameters, randomly extracted from a multidimensional 
sphere, centered on the central value of the parametriza- 
tion, with a radius corresponding to Ict of statistical un- 
certainty. The correlations between the parameters are 
neglected because their total effect, known from Ref. [2J], 
where they have been accounted for in detail, brings a re- 
duction of the final systematic uncertainty because most 
correlations are negative. The dE / dx-'mAxiceA systematic 
uncertainty on each observable is then obtained as the 
standard deviation of the distribution of that observable, 
over the ensemble of likelihood fits performed with dif- 
ferent sets of parameters. This approach is adequate for 
our purposes since the statistical uncertainty is greater 
than or of the same order of the systematic uncertainty. 

The second dominant contribution to the systematic 
uncertainty for B^ — > tt+tt^ comes from the uncertainty 
on the relative efficiency correction, while for i?'^ ^• 
K'^ K~ it comes from the uncertainty in the background 
model, which includes a sizeable component of partially 
reconstructed decays with poorly known branching frac- 
tions. The latter systematic uncertainty is conservatively 
assessed by performing extreme variations of the assumed 
relative contributions of the various modes in the simula- 
tion; the resulting uncertainty is still a factor of 2 lower 
than the uncertainty associated to the dE / dx model. 

Other contributions come from trigger efficiencies, h- 
hadron masses, &-hadron lifetimes and AFs/r^, and 
transverse momentum distribution of the A? baryon. A 
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TABLE II: Measured relative branching fractions of rare modes. Absolute branching fractions were derived by normalizing to 
the current world-average value 5(5" ^ K+tv') = (19.4 ±0.6) x 10~^ and assuming the average values at high energy for the 
production fractions: fa/ fd = 0.282 ± 0.038 [2^. The first quoted uncertainty is statistical; the second is systematic. 

"Mode Relative B Absolute B (10'") Limit (10'") 

B° K+R- I'^BiCy+^-f ^ 0-012 ± 0.005 ± 0.005 0.23 ± 0.10 ± 0.10 [0.05, 0.46] at 90% C.L. 

B°^TT+-K- n^^S^lX"-, ^ 0.008 ± 0.002 ± 0.001 0.57 ± 0.15 ± 0.10 



further systematic uncertainty of the order of 10% is in- 
cluded for the — )■ K^K~ mode to account for a small 
bias of the fitting procedure observed in simulated sam- 
ples. 

The final results are listed in TableHH Absolute branch- 
ing fractions are also quoted, by normalizing to world- 
average values of production fractions and B{B'^ — > 
if+TT") [i^l- The branching fraction measured for the 
B^ — > TT+TT" mode is consistent with and supersedes 
the previous upper limit (< 1.2 x 10^'' at 90% C.L.), 
based on a subsample of the current data It is 

in agreement with predictions obtained with the pQCD 
approach 0, but it is higher than most other the- 
oretical predictions P, d, [25|. The central value for 
B{B^ — ^ K^K^) is the most precise determination of 
this quantity to date, and is in agreement with previ- 
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ous experimental results [Hi, |26| and theoretical predic- 
tions [HQ- It supersedes the previous CDF limit [l3| , 
based on a subsample of the current data. The present 
measurements represent a significant step in reducing a 
source of uncertainty in many theoretical predictions for 
charmless B-decays. The results favor a large annihila- 
tion scenario, which is somewhat unexpected for instance 
in QCDF [27]. 

In summary, we have searched in CDF data for as-yet- 
unmeasured charmless decay modes of neutral h mesons 
into pairs of charged mesons. We report an updated up- 
per limit for the B^ K^K~ mode and the first evi- 
dence for the B^g — >■ tt^tt^ mode and a measurement of 
its branching fraction. 
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